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Abstract. Cosmic neutrino events detected by the IceCube Neutrino Observatory with energy 
> 30 TeV have poor angular resolutions to reveal their origin. Ultrahigh-energy cosmic rays 
(UHECRs), with better angular resolutions at > 60 EeV energies, can be used to check if the 
same astrophysical sources are responsible for producing both neutrinos and UHECRs. We 
test this hypothesis, with statistical methods which emphasize invariant quantities, by using 
data from the Pierre Auger Observatory, Telescope Array and past cosmic-ray experiments. 
We find that the arrival directions of the cosmic neutrinos are correlated with > 100 EeV 
UHECR arrival directions at confidence level ~ 90%. The strength of the correlation decreases 
with decreasing UHECR energy and no correlation exists at energy ~ 60 EeV. A search in 
astrophysical databases within 3° of the arrival directions of UHECRs with energy > 100 EeV, 
that are correlated with the IceCube cosmic neutrinos, resulted in 18 sources from the Surift- 
BAT X-ray catalog with redshift z < 0.06. We also found 3 objects in the Kiihr catalog of radio 
sources using the same criteria. The sources are dominantly Seyfert galaxies with Cygnus A 
being the most prominent member. We calculate the required neutrino and UHECR fluxes to 
produce the observed correlated events, and estimate the corresponding neutrino luminosity 
(25 TeV-2.2 PeV) and cosmic-ray luminosity (500 TeV-180 EeV), assuming the sources are 
the ones we found in the Swift -BAT and Kiihr catalogs. We compare these luminosities with 
the X-ray luminosity of the corresponding sources and discuss possibilities of accelerating 
protons to > 100 EeV and produce neutrinos in these sources. 
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1 Introduction 

The IceCube Neutrino Observatory, the world’s largest neutrino detector, has recently pub¬ 
lished neutrino events collected over 3-year period with energy in the ~ 30 TeV-2 PeV range 
[1]. Shower events, most likely due to v e or v T charge current vN interactions, dominate the 
event list (28 including 3 events with 1-2 PeV energy) while track events, most likely due 
to u /t charge current vN interactions, constitute the rest. Among a total of 37 events about 
15 could be due to atmospheric neutrino (6.61( g) and muon (8.4 ± 4.2) backgrounds. A 
background-only origin of all 37 events has been rejected at 5.7-cr level [1], Therefore a cos¬ 
mic origin of a number of neutrino events is robust. The track events have on average ~ 1° 
angular resolution, but the dominant, shower events have much poorer angular resolution, 
~ 15° on average [1], thus making them unsuitable for astronomy. 

Meanwhile the Pierre Auger Observatory (PAO) [2] and the Telescope Array (TA) [3], 
two of the world’s largest operating cosmic-ray detectors, have recently released UHECR 
data collected over more than 10-year and 5-year periods, respectively. Together they have 
detected 16 events (6 by PAO [2] and 10 by TA [3]) with energies > 100 EeV. The total 
publicly available > 100 EeV events including past experiments is 33. While lower-energy 
cosmic ray arrival directions are scrambled by the Galactic and intergalactic magnetic fields, 
at > 60 EeV energies the arrival directions of UHECRs tend to be much better correlated 
with their source directions and astronomy with charged particles could be realized [4], Few 
degree angular resolution can be achieved at these energies, which is much better than the 
IceCube neutrino shower events and is comparable to the neutrino track events. 

The astrophysical sources of UHECRs with energy > 40 EeV need to be located within 
the so-called GZK volume [5, 6] in order to avoid serious attenuation of flux from them due to 
interactions of UHECRs with photons from cosmic microwave background (CMB) and extra- 
galactic background light (EBL). The astrophysical sources of neutrinos, on the other hand, 
can be located at large distances and still be detected provided their luminosity is sufficiently 
high. However, because of weakly interacting nature of neutrinos and limiting luminosity of 
astrophysical sources, only nearby neutrino sources can be identified, thus making neutrino 
astronomy possible. 
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We explore here a possibility that both UHECRs and IceCube cosmic neutrino events are 
produced by the same astrophysical sources within the GZK volume. Since widely accepted 
Fermi acceleration mechanism of cosmic rays at the sources take place over a large energy 
range, it is natural that the same sources produce > 1 PeV cosmic rays, required to produce 
cosmic neutrinos observed with energies down to ~ 30 TeV, and UHECRs with energy > 40 
EeV. We employ invariant statistical method [7, 8], independent of coordinate systems, in 
order to study angular correlation between cosmic neutrinos and UHECRs. As far we know, 
this is the first attempt to quantify such a correlation between the IceCube neutrino and 
UHECR data sets. Existence of such a correlation can provide clues to the origin of both 
cosmic neutrinos and UHECRs. We search for astrophysical sources within the angular errors 
of UHECRs which are correlated with the neutrino events in order to shed lights on their 
plausible, common origins. Finally we calculate required neutrino and cosmic-ray luminosities 
for the sources to produce observed events, and compare these luminosities with their observed 
X-ray and radio luminosities to check if they are viable sources of both UHECRs and cosmic 
neutrinos. 

The organization of this paper is the following. We describe neutrino and UHECR data 
that we use in section 2 and our statistical method in section 3. The results of our correlation 
study and source search along the directions of the correlated events are given in section 4. 
Section 4 also includes calculation of neutrino and cosmic-ray luminosities of the correlated 
sources from respective fluxes derived using data. We discuss our results and implications of 
our findings in section 5. 


2 IceCube neutrino events and UHECR data 

We consider 35 IceCube neutrino events, collected over 988 days in the ~ 30 TeV-2 PeV range, 
from ref. [1] to study angular correlation with UHECRs. Two track events (event numbers 
28 and 32) are coincident hits in the Ice Top surface array and are almost certainly a pair 
of atmospheric muon background events [1], Therefore we excluded them from our analysis. 
Figure 1 shows sky maps of 35 events in equatorial coordinates with reported angular errors. 
The majority (26) of the events have arrival directions in the southern sky. Among the 9 
northern hemisphere events, only 1 is at a declination Dec > 41° which happens to be a 
shower event [1]. The angular resolutions of the track events are < 1.4° and those of the 
shower events vary between 6.6° and 46.3°. Figure 2 shows sky maps in Galactic coordinates. 

Figures 1 and 2 also show sky maps of available UHECR data with energies > 100 EeV 
(top panel) and > 40 EeV (bottom panel). The PAO and TA collaborations have published 
data with energies above 52 EeV (231 events) [2] and 57 EeV (72 events) [3], respectively. 
Note that the PAO and TA are located in the southern and northern hemisphere, respectively, 
covering a declination range of —90° < Dec < 45° [2] and —10° < Dec < 90° [3]. The angular 
resolutions for the PAO events with energy > 10 EeV is < 0.9° [9] while for the TA events 
with energy > 57 EeV it is between 1.0° and 1.7° [3]. Note that the > 40 EeV data sample 
is incomplete. Only the AGASA experiment has published data above 40 EeV (40 events) 
covering a declination range —10° < Dec < 90° and with angular resolution < 2° [12], Only 
> 100 EeV data are available from the other past experiments: Haverah Park [10, 11], Yakutsk 
[10], Volcano Ranch [10] and Fly’s eye [10]. Note that these were all northern hemisphere 
experiments. We could not include 13 events with energy > 56 EeV from the HiRes experiment 
as the energies of the individual events are not available [13]. 
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Figure 1. Sky maps of the IceCube > 30 TeV cosmic neutrino events with error circles and UHECR 
data in equatorial coordinates. The top panel shows UHECRs with energy > 100 EeV and the bottom 
panel shows all available data with energy > 40 EeV. The black dotted line is the Galactic plane. 


We list in Table 1 all available UHECR data with energy > 100 EeV. This includes 
6 events from Haverah Park, 1 event from Yakutsk, 8 events from AGASA, 1 event from 
Volcano Ranch, the highest-energy (320 EeV) event from Fly’s eye, 6 events from PAO and 
10 events from TA. We use this list and sublists with PAO and TA data separately to study 
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Figure 2. The same as Figure 1 but in Galactic coordinates. 


correlation with cosmic neutrino events. In addition to > 100 EeV data, we also explore 
energy-dependence of correlation by choosing different energy cuts, > 80 EeV and > 60 EeV, 
in UHECR data from the PAO, TA and AGASA . Above 80 EeV (60 EeV) there are 22 (136) 
UHECR events from PAO, 20 (60) from TA and 11 (22) from AGASA. We note that energy 
calibration across the experiments can vary by as much as ~ 30% (see, e.g., ref. [14]). We 
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Experiment 

Reference 

Energy (EeV) 

RA (°) 

Dec (°) 

Haverah Park 

[10] 

101 

201 

71 

Haverah Park 

[10] 

116 

353 

19 

Haverah Park 

[10] 

126 

179 

27 

Haverah Park 

[10] 

159 

199 

44 

Haverah Park 

[11] 

123 

318.3 

3.0 

Haverah Park 

[11] 

115 

86.7 

31.7 

Yakutsk 

[10] 

110 

75.2 

45.5 

AGASA 

[12] 

101 

124.25 

16.8 

AGASA 

[12] 

213 

18.75 

21.1 

AGASA 

[12] 

106 

281.25 

48.3 

AGASA 

[12] 

144 

241.5 

23.0 

AGASA 

[12] 

105 

298.5 

18.7 

AGASA 

[12] 

150 

294.5 

-5.8 

AGASA 

[12] 

120 

349.0 

12.3 

AGASA 

[12] 

104 

345.75 

33.9 

Volcano Ranch 

[10] 

135 

306.7 

46.8 

Fly’s eye 

[10] 

320 

85.2 ±0.5 

48.0±i;| 

Pierre Auger 

[2] 

108.2 

45.6 

-1.7 

Pierre Auger 

[2] 

127.1 

192.8 

-21.2 

Pierre Auger 

[2] 

111.8 

352.6 

-20.8 

Pierre Auger 

[2] 

118.3 

287.7 

1.5 

Pierre Auger 

[2] 

100.1 

150.1 

-10.3 

Pierre Auger 

[2] 

118.3 

340.6 

12.0 

Telescope Array 

[3] 

101.4 

285.74 

-1.69 

Telescope Array 

[3] 

120.3 

285.46 

33.62 

Telescope Array 

[3] 

139.0 

152.27 

11.10 

Telescope Array 

[3] 

122.2 

347.73 

39.46 

Telescope Array 

[3] 

154.3 

239.85 

-0.41 

Telescope Array 

[3] 

162.2 

205.08 

20.05 

Telescope Array 

[3] 

124.8 

295.61 

43.53 

Telescope Array 

[3] 

135.5 

288.30 

0.34 

Telescope Array 

[3] 

101.0 

219.66 

38.46 

Telescope Array 

[3] 

106.8 

37.59 

13.89 


Table 1 . Available UHECR data with energy > 100 EeV from various experiments, 
discuss this issue in Sec. 4.2. 

3 Statistical method and data analyses 

To study correlation between cosmic neutrinos and UHECRs, we map the Right Ascension 
and Declination (RA, Dec) of the event directions into unit vectors on a sphere as 

x = (sin 6 cos 0, sin 6 sin 4>, cos 0 ) T , 

where (f> = RA and 9 = tt/ 2 — Dec. Scalar product of the neutrino and UHECR vectors 
(^neutrino ' ®uhecr) therefore is independent of the coordinate system. The angle between the 
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two vectors 


7 — COS (^neutrino ' ^UHECr)j (3.1) 

is an invariant measure of the angular correlation between the neutrino and UHECR arrival 
directions [7, 8 ]. Following ref. [7] we use a statistic made from invariant 7 for each neutrino 
direction Xi and UHECR direction x 3 pair as 

5x 2 i = min, ( 7 ^ /Sr/f ), (3.2) 

which is minimized for all j. Here <$ 7 * is the 1-a angular resolution of the neutrino events. 
We use the exact resolutions reported by the IceCube collaboration for each event [1]. 

A value Sxf < 1 is considered a “good match” between the i-th neutrino and an UHECR 
arrival directions. We exploit distributions of all S'x'j statistic to study angular correlation 
between IceCube neutrino events and UHECR data. The distribution with observed data 
giving a number of “hits” or Wits with <5y 2 < 1 therefore forms a basis to claim correlation. 
Note that in case more than one UHECR directions are within the error circle of a neutrino 
event, the <5y 2 value for UHECR closest to the neutrino direction is chosen in this method. 

We estimate the significance of any correlation in data by comparing Wits with corre¬ 
sponding number from null distributions. We construct two null distributions, in one case 
we randomize only the RA of UHECRs, keeping their Dec the same as in data; and in the 
second case we also randomize Dec according to the zenith-angle depended sky exposure of 
the UHECR experiments [36], affecting the declination distributions of UHECR data. We call 
these two null distributions as the semi-isotropic null and exposure-corrected null, respectively. 
The semi-isotropic null is a quick-way to check significance while the exposure-corrected null is 
accurate when information on particular experiments are available. In both cases we perform 
100,000 realizations of drawing random numbers to assign new RA and Dec values for each 
event to construct dy 2 distributions in the same way as done with real data. We find that 
the two null distributions are in good agreement with each other in most cases. 

We calculate statistical significance of correlation in real data or p-value (chance proba¬ 
bility) using frequentists’ approach. We count the number of times we get a random data set 
that gives equal or more hits than the Wits hi real data within <5y 2 < 1 bin. Dividing this 
number with the total number of random data sets generated ( 100 , 000 ) gives us the p-value. 
We cross-check this p-value by calculating the Poisson probability of obtaining Wits within 
dx 2 < 1 bin given the corresponding average hits expected from the null distribution. These 
two chance probabilities are in good agreement. 

4 Results 

4.1 Correlations between neutrinos and UHECRs 

We apply our statistical method separately to the PAO and TA data, to a combination of 
the both and to all available UHECR data above 100 EeV from all experiments (see Table 
1). The results are shown in the histograms of Figure 3. The counts in the 0 < h \ 2 < 1 bins 
for the blue, filled histograms correspond to the number of correlated neutrino events with 
UHECRs. Counts in other bins are due to distant pairs of the neutrino events and UHECRs 
and are uninteresting for us. The counts for the red (green), open histogram in the same 
bins correspond to the expected number of correlated neutrino events from the semi-isotropic 
null ( exposure-corrected null), after averaging over 100,000 simulated data sets with random 
UHECR positions. Both null distributions give similar results. 
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PAO > 100 EeV 


TA > 100 EeV 




TA+PAO > 100 EeV 


All > 100 EeV 




5 x 2 5 '/- 2 


Figure 3. Distributions of <5\ 2 found in observed data (blue, filled histograms) and in simulated data 
corresponding to the semi-isotropic null (red, open histograms) and the exposure-corrected null (green, 
open histograms). The histograms have been truncated at S\ 2 = 10 for better display. Significances 
(p-values) have been calculated for the 0 < Sx 2 < 1 bins. 


Figure 3 also shows p -values or the probability of finding the correlated events (0 < 
<5 \' 2 < 1) in observed data as a fluctuation of the randomly distributed UHECRs in the sky. 
The probability 1 — p « 90% is the confidence level (CL) that the IceCube neutrino events 
and all available UHECR data with energy > 100 EeV are correlated. A correlation with 
similar CL exists between the neutrino and PAO-only data sets and between the neutrino 
and TA-only data sets. The Poisson probability of obtaining Ah; ts = 7 in PAO data when 
4 are expected from the semi-isotropic null distribution and ; ts = 6 in TA data with 3.8 
expected from the same null distribution are ~ 0.06, in very good agreement with our p- values. 
Similarly, for the combined data set of all UHECRs > 100 EeV, Vhits = 12 when expected 
value is 8.8 corresponds to a Poisson probability of ~ 0.07, again in very good agreement with 
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Figure 4. Comparisons between the Ahits distributions in the 6\ 2 < 1 bins of Fig. 3 obtained 
from the semi-isotropic (red solid lines) and exposure-corrected (green solid lines) null distributions. 
Also plotted are the Poisson distributions (dotted lines) for the average values of the respective null 
distributions in the Sx 2 < 1 bins of Fig. 3. The vertical lines are the observed IVhits values in data. 


our p-value. 

We remind the readers that the counts in the Sx' 2 distributions with TA+PAO data is 
not the algebraic sum of the counts in the distributions with TA and PAO data separately. 
This is because our Sx 2 statistic choose the nearest UHECR data point even if more than one 
are present within the error circle of a neutrino event. The same is true for distribution with 
all UHECRs. We list the correlated events in Table 2 against the neutrino event numbers in 
ref. [1], Note that we list all UHECRs giving Sx 2 < 1 in the table. There are 7 UHECRs 
which are correlated with 2 or more neutrino events. None of the correlated neutrinos are 
PeV neutrino events. 

Figure 4 shows a comparison between the two null distributions for UHECRs with 










v event no. [1] 

$x 2 

Energy (EeV) 

RA (°) 

Dec (°) 

Experiment 

1 

0.41 

108.2 

45.6 

-1.7 

PAO 


0.95 

106.8 

37.59 

13.9 

TA 

2 

0.97 

150 

294.5 

-5.8 

AGASA 

11 

0.10 

100.1 

150.1 

-10.3 

PAO 

16 

0.006 

127.1 

192.8 

-21.2 

PAO 

17 

0.77 

144 

241.5 

23.0 

AGASA 

21 

0.55 

111.8 

352.6 

-20.8 

PAO 

24 

0.78 

101.4 

285.74 

-1.7 

TA 


0.97 

150 

294.5 

-5.8 

AGASA 

25 

0.06 

150 

294.5 

-5.8 

AGASA 


0.07 

101.4 

285.74 

-1.7 

TA 


0.10 

135.5 

288.3 

0.34 

TA 


0.12 

118.3 

287.7 

1.5 

PAO 


0.58 

105 

298.5 

18.7 

AGASA 


0.62 

123 

318.3 

3 

Haverah Park 

29 

0.18 

124.8 

295.6 

43.52 

TA 

31 

0.35 

101 

201 

71 

Haverah Park 

33 

0.34 

118.3 

287.7 

1.5 

PAO 


0.40 

135.5 

288.3 

0.34 

TA 


0.74 

101.4 

285.74 

-1.7 

TA 


0.84 

105 

298.5 

18.7 

AGASA 

34 

0.20 

104 

345.75 

34 

AGASA 


0.22 

135 

306.7 

46.8 

Volcano Ranch 


0.25 

122.2 

347.7 

39.46 

TA 


0.34 

118 

340.6 

12 

PAO 


0.34 

124.8 

295.61 

43.53 

TA 


0.36 

105 

298.5 

18.7 

AGASA 


0.45 

123 

318.3 

3 

Haverah Park 


0.47 

116 

353 

19 

Haverah Park 


0.50 

120 

349 

12.3 

AGASA 


0.55 

120.3 

285.5 

33.62 

TA 


0.71 

134 

281.25 

48.3 

AGASA 


Table 2. IceCube cosmic neutrino events correlated with UHECRs above 100 EeV. 


energy > 100 EeV in Fig. 3 by using the A^its within the <5y 2 < 1 bin from simulations for 
both the null distributions. The two null distributions agree well in all cases except for the 
combined analysis of the TA and PAO data. A comparison with Poisson distribution with 
frequency for the corresponding cases are also shown. For PAO UHECRs >100 EeV the two 
null distributions follow the respective Poisson distributions but not for the other cases. The 
black vertical line represents the observed N^ ts . 

We do the same statistical analysis with UHECRs above 80 EeV. The results are shown 
in Figure 5. There are 6 new correlations in the PAO-only data. Only 10 is expected from 
our null distributions as compared to 13 total in data. This reduces correlation between the 
IceCube and PAO data to ~ 84% — 88% CL. A list of UHECRs from the PAO correlated with 
the neutrino events is given in Table 3. Note that the 2 PeV neutrino event (event number 
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PAO > 80 EeV 


TA > 80 EeV 




TA+PAO > 80 EeV 


TA+PAO+AGASA > 80 EeV 




5x 2 5 * 2 


Figure 5. The same as in Figure 3 but for UHECRs with energy > 80 EeV. 


35), the highest observed energy so far [1], is now correlated with a 89.1 EeV PAO event. The 
other two neutrino events with energy ~ 1 PeV each still remain uncorrelated with UHECRs 
with energy > 80 EeV. There are 6 UHECRs in Table 3 which are correlated with more than 
one neutrino events. In particular an 82 EeV PAO event is correlated with 4 neutrino events. 

Lowering the UHECR energy lower limit to 80 EeV adds 2 new correlated events in 
the case of TA-only data with a total of 8 as compared to ~ 5.6 expected from both the null 
distributions, giving a CL of ~ 90%. Combining the PAO and TA data results in similar 
significance as obtained from individual data sets. Combining the PAO, TA and AGASA 
data reduces the significance of correlation. 

Further lowering the UHECR energy lower limit to 60 EeV gives no significant correlation 
between the IceCube cosmic neutrino data and UHECR data. Figure 6 shows that the number 
of correlated events in data is very similar to those expected from the null distributions in all 
cases. Such a loss of significance is expected when there is no real correlation between the 
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v event no. [1] 


Energy (EeV) 

RA (°) 

Dec (°) 

1 

0.41 

108.2 

45.6 

-1.7 

2 

0.002 

80.9 

283.7 

-28.6 

7 

0.85 

83.8 

26.8 

-24.8 

11 

0.10 

100.1 

150.1 

-10.3 

15 

0.6 

82.3 

287.7 

-64.9 

16 

0.006 

127.1 

192.8 

-21.2 


0.52 

84.7 

199.7 

-34.9 

21 

0.55 

111.8 

352.6 

-20.8 

21 

0.6 

83.8 

26.8 

-24.8 

22 

0.85 

80.9 

283.7 

-28.6 

24 

0.77 

80.9 

283.7 

-28.6 

25 

0.095 

80.2 

283.7 

-28.6 


0.12 

118.3 

287.7 

1.5 


0.61 

82 

299 

19.4 


0.62 

80.2 

271.1 

19.0 


0.67 

81.4 

308.8 

16.1 

33 

0.34 

82 

287.7 

1.5 


0.95 

82 

299 

19.4 

34 

0.22 

81.4 

308.8 

16.1 


0.34 

82 

299. 

19.4 


0.34 

118.3 

340.6 

12 


0.6 

89 

349.9 

9.3 

35 

0.96 

89.1 

218.8 

-70.8 


Table 3. 


IceCube cosmic neutrino events correlated with UHECRs detected by PAO above 80 EeV. 


data sets. 

4.2 Energy calibrated events 

As we noted earlier, energy calibration among different UHECR experiments is a widely- 
debated issue. If UHECR flux is uniform over the whole sky and each experiment measures 
the same primary particles’ energy then the number of UHECR events should be proportional 
to the exposures of the experiments. Therefore it is difficult, in particular, to reconcile the 
10 TA events at > 100 EeV compared to the 6 from PAO, which has 20 times more exposure 
than TA. There are many energy rescaling procedure suggested among experiments (see, 
e.g., refs. [15-17]) to bring their respective measured fluxes close to each other, mostly at 
the “ankle” regime. Even these procedures cannot reconcile number of events, after exposure 
corrections, above 100 EeV among different experiments. It is plausible that the energy 
rescaling factors themselves are energy dependent, differing from the ankle regime to the 
GZK regime. Reconstructing such energy-dependent rescaling factors is beyond the scope of 
this paper. We hope the experimental collaborations will provide such factors in future. 

To illustrate the energy rescaling effect on our correlation study, we adopt a recent 
procedure in ref. [15] which is based on a joint PAO and TA analysis of UHECRs from an 
overlapping region in the sky. We decrease the energies of the TA events by 25% but keep 
the energies of the PAO events unchanged [15]. So the number of UHECRs events above 
100 EeV from TA is now 4, out of which only 2 correlates with the neutrino events in the 
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Figure 6. The same as in Figure 3 but for UHECRs with energy > 60 EeV. 


<5y 2 < 1 bin (see Fig. 7). Interestingly, however, there are now 6 counts in the 1 < <5y 2 < 2 
bin which corresponds to a Poisson probability of 0.0108 according to the semi-isotropic null. 
A combined analysis of the TA and PAO data considering the above energy rescaling, gives 
no significant correlation with neutrino data. 


4.3 Astrophysical source search 

We search for astrophysical source candidates for UHECRs which are correlated with IceCube 
cosmic neutrino events, assuming both are produced by the same sources. We use data from 
Tables 2 and 3 for this purpose. The experimental angular resolution of the UHECRs is of the 
order of 1°. However, Galactic and intergalactic magnetic held can deflect them by more than 
a few degrees from their source directions. The deflection angle in the intergalactic random 
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TA > 100 EeV 


TA+PAO > 100 EeV 




Figure 7. The same as in Figure 3 upper-right and bottom-left panels but the energies for the TA 
events have been reduced by 25% compared to the PAO events. 


56i G « 1.1°Z 


(4.1) 


magnetic field [18] with strength B r d m and coherence length A co h is 

_1 ( B ^\ ( D V /2 ( A coh y /2 
100 EeV ) V lnG /V 200 M P c / VlOOkpcy 

where Z and E cr are the charge and energy of the UHECR and D is the distance to the 
source. The deflection angle in a small-scale Galactic random magnetic field, using Eq. (4.1) 
with -B r dm = lyuG, A co h = 100 pc is much smaller, 59q ~ 0 . 2 °Z, for E cr = 100 EeV and 
D = 10 kpc. However, the deflection angle in the large-scale regular component of the 
Galactic magnetic field in the disk and in the halo can be larger, ~ l°-3° [19, 20]. Hereafter 
we assume that UHECRs with energy > 80 EeV are dominantly protons 1 . We also chose a 
conservative source search region of 3° around the directions of UHECRs which are correlated 
with cosmic neutrino events. 

We also limit our source search within a comoving volume with its radius set by the GZK 
effect of Puhecr + 7CMB interactions and corresponding energy losses by UHECR protons 
[5, 6 ]. A crude estimate of the mean-free-path for this interaction can be obtained from the 
number density of CMB photons with 2.73 K temperature in the local universe, which is 

1.32 x 10 4 (e/meV ) 3 . _ 3 


en(e) = 


cm 


(4.2) 


exp [4.25(e/meV)[ — 1 

Thus the number density is 244 cm” 3 at the peak photon energy e = 2.82 £;b( 2.73 K) = 
0.66 meV, where ks = 8.62 X 10” 5 eV K” 1 is the Boltzmann constant. A parametrization 
of the UHECR proton’s mean-free-path, using delta function approximation of the py cross 
section, is given by 

-3 


A„ ss 245.76 


E p 

100 EeV 


exp 


0.42 


E p 

100 EeV 


- 1 


Mpc, 


(4.3) 


1 Note that the mass composition measurement by the PAO collaboration by using shower maxima, which 
favors heavy nuclei as primaries, is done up to an energy ~ 60 EeV only [22]. 
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which reproduces results from numerical calculations with accurate treatment [ 21 ] within 
~ 10% in the ~ 60-200 EeV range. For reference, \ p = 539, 247, 138 and 26 Mpc at 
E p = 60, 80, 100 and 200 EeV, respectively. We search for astrophysical sources within 
redshift z = 0.06, which corresponds to a luminosity distance df = 270.4 Mpc in ACDM 
cosmology with Ho = 69.6 km s -1 Mpc -1 , Qm = 0.286 and = 0.714 [23]. The proper 
distance, d p = di,/( 1 + z) 2 = 241 Mpc, is similar to \ p at 80 EeV. 

We have used the Swift- BAT 70 month X-ray source catalog [24] to search for astrophys¬ 
ical sources which are correlated with UHECR and cosmic neutrino events. In 70 months 
of observations, the catalog includes 1210 objects of which 503 objects are within redshift 
< 0.06. Out of these 503 X-ray selected objects at least 18 are simultaneously correlated with 
the neutrino events and UHECRs above 100 EeV, see Table 4. The correlated X-ray sources 
are all Seyfert (Sy) galaxies except ABELL 2319 which is a galaxy cluster (GC). The X-ray 
luminosity of these sources vary between L \ ~ 10 43 -10 45 erg s -1 , with Cygnus A the most 
luminous of all. Note that the PAO collaboration has also found an anisotropy at ~ 98.6% CL 
in UHECRs with energy > 58 EeV and within ~ 18° circles around the AGNs in Swift- BAT 
catalog [24] at distance < 130 Mpc and X-ray luminosity Lx > 10 ' 44 erg s -1 [2], Our list in 
Table 4 includes NGC 1142 which is also one of the five sources that dominantly contribute 
to the anisotropy found in the PAO data [ 2 ], 

In another correlated source search we have used bright extragalactic radio sources with 
flux density > 1 Jy at 5 MHz from the Kiihr catalog [25]. It has 61 sources within known 
redshift < 0.06. Only 3 sources from this catalog are correlated simultaneously with IceCube 
neutrinos and UHECRs above 100 EeV, see Table 4. Two of these sources are Seyfert galaxies 
and the third one is a galaxy cluster. There are two common sources, that are correlated with 
both neutrinos and UHECRs, between the Swift-BAT and Kiihr catalogs. These are NGC 
1068 and PKS 2331-240. Both of them are Seyfert galaxies. 

We have also searched the first AGN catalog (1LAC) published by the Fermi- LAT col¬ 
laboration [26] for possible correlations with neutrino and UHECR arrival directions but did 
not find any. 

It is interesting note that the cosmic neutrino events (nos. 2, 12, 14, 15 and 36) which 
are strongly correlated with the Fermi Bubbles 2 [27, 28], except for event no. 2, do not appear 
in Table 4. This could be a hint to possible extragalactic [29] and Galactic [30] components 
in the neutrino event data. 


4.4 Neutrino and UHECR luminosities for correlated events 


After searching for astrophysical sources correlated with both IceCube cosmic neutrino events 
and UHECRs, we calculate their corresponding fluxes required to produce observed events. 
First, we describe our point-source neutrino flux calculation method. We assume a power-law 
flux of the following form 

J Va {Ev) = A u a ( 100 ^ eV ) , (4.4) 

which is the same for all 3 flavors: a = e, //, r. We estimate the normalization factor from 
the number of neutrinos events N v of any flavor 3 as 


^ Je vI d-Ej, A e ff a (Ev) ( ] 0 (VTeV) 

2 The centers of the error circles within the Fermi bubbles’ contours 

3 Here we tacitly assume that flavor identification is not efficient. 


(4.5) 
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Neutrino 
Event # 


UHECR 

Swift X-ray Source Catalog [24] 


RA 

Dec 

Experiment 

Name 

z 

Type 

1 

45.6 

-1.7 

PAO 

NGC 1142 

0.0289 

Sy2 





NGC 1194 

0.0136 

Syl 





MCG +00-09-042 

0.0238 

Sy2 





NGC 1068 

0.0038 

Sy2 

11 

150.1 

-10.3 

PAO 

2MASX J10084862-0954510 

0.0573 

Syl .8 

17 

241.5 

23 

AGASA 

2MASX J16311554+2352577 

0.0590 

Sy2 

29, 34 

295.6 

43.52 

TA 

2MASX J19471938+4449425 

0.0539 

Sy2 





ABELL 2319 

0.0557 

GC 





Cygnus A 

0.0561 

Sy2 

21 

352.6 

- 20.2 

PAO 

PKS 2331-240 

0.0477 

Sy2 

2, 24, 25 

294.5 

-5.8 

AGASA 

2MASX J19373299-0613046 

0.0103 

Syl.5 

34 

340.6 

12 

PAO 

MCG +01-57-016 

0.0250 

Syl .8 





MCG +02-57-002 

0.0290 

Syl.5 





UGC 12237 

0.0283 

Sy2 


349.0 

12.3 

AGASA 

NGC 7479 

0.0079 Sy2/Liner 





2MASX J23272195+1524375 

0.0457 

Syl 





NGC 7469 

0.0163 

Syl.2 


352.6 

- 20.2 

Haverah Park 

NGC 7679 

0.0171 

Sy2 

Neutrino 
Event # 


UHECR 

Ruhr Radio Source Catalog [25] 


RA 

Dec 

Experiment 

Name 

z 

Type 

1 

45.6 

-1.7 

PAO 

NGC 1068 

0.0038 

Sy2 

21 

352.6 

- 20.8 

PAO 

PKS 2331-240 

0.0477 

Sy2 

34 

340.6 

12 

PAO 

NGC 7385 

0.0255 

GC 


Table 4. Sources correlated with UHECRs and neutrino events simultaneously. 


where T is IceCube lifetime and A e g >a is effective area for different flavors. We use T = 988 
days for IceCube 3-year data release [1] as in our correlation analysis and the following 
parametrization, correct within ~ 10% uncertainty, of the effective areas [31] 


AeS,e 


Aeff ,/i 


A 


eff ,r 


[1.26 x 1(T 5 (E v / TeV) 2 - 64 - 0.017] nr 2 , 
[0.459 (E v / TeV) 0 ' 5 - 1.109] nr 2 , 

[3.6 x 1(T 6 (E^/TeV) 2 ' 64 - 0.0142] m 2 , 
[0.389 (E u / TeV) 0 ' 5 - 1.868] nr 2 , 

[7.267 x 1(T 6 (Eu/ TeV) 2 - 64 - 0.0175] m 2 , 
[0.5069 (E v / TeV) 0 ' 5 - 3.092] nr 2 , 


25 TeV < E u < 100 TeV 
E u > 100 TeV 

25 TeV <E U < 100 TeV 
Eu > 100 TeV 

25 TeV <E V < 100 TeV 
E v > 100 TeV. 
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For the limits of the integral in Eq. (4.5) we set E v \ = 25 TeV and E u 2 = 2.2 PeV, 
reflecting uncertainty in energy estimate reported by the IceCube collaboration [1]. 

We use neutrino flux to calculate neutrino luminosity of the corresponding source as 

fEv2(l+z) 

L v = 4vt4 V / dE v E v J Vja (E v ). (4.7) 

q, J E v i(l +z) 


These luminosities are listed in Table 5. We use two values of k. the choice k = 2.3 is motivated 
by fit to IceCube data assuming an isotropic distribution of events [1, 32] and the choice k = 
2.1 is motivated by the cosmic-ray spectrum expected from Fermi acceleration mechanisms. 
Note that N v = 3 for 2MASX J19373299-0613046, N„ = 2 for 2MASX J19471938+4449425 
and N u = 1 for all other sources. Neutrino luminosities listed in Table 5 are within a factor 
5 of the corresponding X-ray luminosities of the sources. For the radio sources, neutrino 
luminosity far exceeds the corresponding radio luminosity. 

We also calculate UHECR flux from the observed events by using exposure for the 
respective experiments. We use a power-law form for the observed UHECR flux above 28.8 
EeV break energy as [33] 

E \ ~ 4 ' 3 

; E cr > 28.8 EeV, (4.8) 


'4uhecr(7?cr) — ^4uhecr 


and derive the normalization factor as 


A 


uhecr 


Ai sheer 


Su(S) rE CT 2 

fi J-Ecrl 


( EeV ) 


-4.3 ' 


(4.9) 


Here S is the total integrated exposure, as mentioned in ref. [34], H is the solid angle of the de¬ 
tector and oj(6) is the relative exposure for particular declination angle <5. For reference, we use 
for PAO, SpAo = 66, 000 km 2 yr sr and Hpao = 1-657T sr [2]; for TA, SpA = 3, 690 km 2 yr sr 
and Hta = 0.85-7T sr [3]; for AG AS A, Hagasa = 1,000 km 2 yr sr and Hagasa = 0 .597T sr 
[35]. We do not use the Haver ah Park event that is correlated with a neutrino event. We 
calculate co(6) from ref. [36] but adapt it for different experiments by using their respective 
geographical locations and zenith angle ranges. For the lower and upper limits of integration 
in Eq. (4.9), we use E cr \ = 80 EeV and E CI 2 = 180 EeV, allowing a 20% uncertainty for the 
100 EeV threshold energy used to search correlation and 150 EeV maximum energy found for 
a UHECR correlated with a neutrino event and a source (see Table 2). 

In order to calculate cosmic-ray 4 luminosity of the sources, first we note that if cosmic 
neutrinos detected by IceCube are coming from the same source that are correlated with 
UHECRs, then the cosmic-ray flux in Eq. (4.8) needs to be extrapolated down to ~ 500 TeV 
which is required to produce ~ 30 TeV neutrinos. Second, the flux in Eq. (4.8) needs to be 
corrected for GZK suppression above ~ 40 EeV. Given a cosmic-ray proton luminosity L cr 
between the generation energies E' crl = 500 TeV and E' cr2 = 180 EeV with oc E'~ K spectrum 
in a source at redshift z , the cosmic-ray flux on the Earth is [37, 38] 


Jcr(E c 


A cr (l + z) (k- 2)(E' cvl E> v2 r- 2 (dE’\ 

4v Td\ E' c r ;- 2 - E'£- 2 cr \dE c J ' 


(4.10) 


The cosmic-ray energy at the source and on the Earth are related through various energy 
losses [37]. Following ref. [38] we have plotted cosmic-ray flux in Fig. 8 using Eq. (4.10) for 

4 We assume they are dominantly protons. 
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Source name 

L x (10 44 erg/s) 

L v (10 44 

erg/s) 

L cr (10 44 

erg/s) 


/L r (10 41 erg/s) 

K = 2.1 

= 2.3 

K = 2.1 

= 2.3 

NGC 1142 

1.58/0.012(74 GHz) 

0.95 

1.0 

0.7 

5.4 

NGC 1194 

0.12/0.00012(1.4 GHz) 

0.2 

0.2 

0.04 

0.2 

MCG +00-09-042 

0.17/0.0043(1.4 GHz) 

0.64 

0.71 

0.3 

2.1 

NGC 1068 

0.031/0.0034(31.4 GHz) 

0.016 

0.017 

0.001 

0.007 

2MASX J10084862-0954510 

1.04/0.0028(1.4 GHz) 

3.9 

4.32 

44 

578 

2MASX J16311554+2352577 

0.79/0.0048(1.4 GHz) 

4.1 

4.6 

1600 

22000 

2MASX J19471938+4449425 

1.66/0.0045(1.4 GHz) 

6.8 

7.6 

211 

26000 

ABELL 2319 

1.78/0.0046(1.4 GHz) 

3.7 

4.1 

270 

3500 

Cygnus A 

11.2/314(14.7 GHz) 

3.7 

4.1 

290 

3700 

PICS 2331-240 

0.81/1.32(31.4 GHz) 

2.6 

2.9 

9.5 

102 

2MASX J19373299-0613046 

0.055/0.0012(1.4 GHz) 

0.24 

0.26 

1.3 

7.3 

MCG +01-57-016 

0.23/0.0026(1.4 GHz) 

0.71 

0.78 

0.5 

3.6 

MCG +02-57-002 

0.25/0.00084(1.4 GHz) 

0.95 

1.1 

1.0 

7.5 

UGC 12237 

0.23/0.0011(1.4 GHz) 

0.91 

1. 

0.9 

6.6 

NGC 7479 

0.029/0.04(22 GHz) 

0.07 

0.08 

0.3 

1.4 

2MASX J23272195+1524375 

0.51/0.24(1.4 GHz) 

2.4 

2.7 

280 

2900 

NGC 7469 

0.4/0.0056(365 MHz) 

0.3 

0.3 

2.2 

14 

NGC 7679 

0.1/0.00033(1.4 GHz) 

- 

- 

- 

- 

NGC 1068 

0.031/0.0034(31.4 GHz) 

0.016 

0.017 

0.001 

0.007 

PICS 2331-240 

0.81/1.32(31.4 GHz) 

2.6 

2.9 

9.5 

102 

NGC 7385 

- /0.17(31.4 GHz) 

0.7 

0.8 

0.5 

4.0 


Table 5. Neutrino (25 TeV-2.2 PeV) and cosmic-ray (500 TeV-180 EeV) luminosities required for the 
correlated sources in Table 4 to produce observed data. Also listed are Swift-BAT X-ray luminosity 
[24] radio luminosity for these sources, with corresponding radio frequencies in parentheses. 


L cr = 10 44 erg s _1 and for various redshift in the range 0.01 < z < 0.06. We have also used 
two different values for k as we did for neutrino flux calculation. 

Figure 8 provides a map to estimate cosmic-ray luminosity of the sources listed in Table 
4 which are correlated with UHECR and neutrino events. The UHECR flux in Eq. (4.8), 
calculated from data, corresponds to a point in Fig. 8 at E CI ~ 80 EeV. We estimate the 
source luminosity L cr by equating this flux to the expected flux in Eq. (4.10) at 80 EeV for 
the redshift of a given source. These luminosities are listed in Table 5. Note that except for 
NGC 1068 (z = 0.0038), NGC 1194 (z = 0.0136) and MCG +00-09-042 (z = 0.0238) the 
cosmic-ray luminosity with n = 2.1 is comparable or higher than the neutrino luminosity for 
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Ecr(eV) 


Figure 8. Expected UHECR flux on the Earth from sources at different redshift 0.01 < z < 0.06 
but with fixed luminosity L p = 10 44 erg s _1 in the 500 TeV to 200 EeV range. 

all sources. The cosmic-ray luminosity exceeds the X-ray or radio luminosities for all sources 
except for NGC 1142 and NGC 1194, in case k = 2.1. 


5 Discussion and Outlook 

We have investigated whether the arrival directions of cosmic neutrinos, detected by IceCube 
[1], with energy ~ 30 TeV-2 PeV are correlated with the arrival directions of UHECRs with 
energy > 100 EeV. In order to test correlation we have used an invariant statistic, called the 
minimum 5x 2 [7], which is constructed from the angle between two unit vectors corresponding 
to the directions of the neutrino events and UHECRs, and weighted by the angular resolutions 
of the neutrino events. We have evaluated the significance of any correlation by using Monte 
Carlo simulations of randomly generated UHECR directions and comparing with data. We 
found that IceCube cosmic neutrinos are correlated with UHECRs with energy > 100 EeV 
with significance at 90% CL. The significance, however, decreases with decreasing energy of 
UHECRs, leaving no correlation at an energy threshold of 60 EeV. To take into account 
trial factor, since we searched for correlation with Aerial = 3 UHECR energy thresholds, we 
calculate post-trial p-value as p pos t-trial = 1 — (1 — p B i g n ai ) 1/ ' jVtriaJ = 0.27, with p s i gna l = 0.1 
that we found in data. 

We have searched for astrophysical sources in the Swift -BAT X-ray catalog [24], the 
Kiihr radio source catalog [25] and Fermi-LAT 1LAC AGN catalog [26] within 3° error circles 
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of the > 100 EeV UHECRs which are correlated with cosmic neutrino events, assuming the 
UHECRs are protons. We made a cut in redshift, z < 0.06, while searching for sources in the 
catalogs. This corresponds to a proper distance of 241 Mpc, similar to the mean-free-path 
of an 80 EeV proton in the CMB. The choice of 3° error circle is motivated by deflection of 
UHECR protons in the intergalactic and Galactic magnetic fields. We found that 18 sources 
from the Swift -BAT X-ray catalog and 3 sources from the Kiihr radio source catalog are 
within 3° error circles of the UHECRs that are correlated with cosmic neutrinos. Except for 
ABELL 2319 and NGC 7385 which are galaxy clusters, the rest of the sources are Seyfert 
galaxies with Cygnus A being the most well known. Our finding is consistent with that of 
the PAO collaboration who found significant correlation between UHECR arrival directions 
and Seyfert galaxies in the Swift -BAT X-ray catalog [2], We did not find any source from the 
Fermi -LAT 1LAC AGN catalog fitting our search criteria. 

Estimates of the neutrino and UHECR fluxes for the correlated events were used to 
calculate corresponding 25 TeV-2.2 PeV neutrino luminosity and 500 TeV-180 EeV cosmic- 
ray luminosity under the hypothesis that both originated from the sources we found in the 
Swift -BAT and Kiihr catalogs. The neutrino luminosities are of the same order as the X-ray 
luminosities of the sources. The cosmic-ray luminosities, depending on the source spectrum, 
are comparable or higher than both the neutrino and X-ray luminosities. Comparison between 
the nonthermal X-ray luminosity with the cosmic-ray or neutrino luminosity gives a possibility 
that the energy in X-ray producing electrons can be compared to that of cosmic-ray protons, 
both accelerated at the sources. 

Acceleration of UHECRs near the central black holes of AGNs was proposed over 20 years 
ago [39, 40]. Interactions of these UHECRs with UV and X-ray photons could produce high- 
energy neutrinos [39, 41]. Seyfert galaxies are radio-quite AGNs and do not have strong jets, 
although parsec scale jets in them have been observed in the last decade [42-44], Collisions 
between blobs in this jet and formation of shocks may lead to acceleration of protons to 
an energy at least up to 10 18 eV, with subsequent photomeson interactions producing high- 
energy neutrinos [45]. Acceleration of heavy nuclei and subsequent gamma-ray and neutrino 
production in radio-quite AGNs have also been discussed [46]. Predictions have also been 
made for GeV-TeV gamma-ray emission from UHECR interactions in Cygnus A [47], which 
is also a powerful radio galaxy (3C 405). 

If a fraction ijx of the X-ray luminosity, Lx = 10 44 L 44 erg s _1 , of the Seyfert galaxies is 
nonthermal then one can estimate the energy density in magnetic field in the X-ray emitting 
region as B 2 /8n = rjxLx/^ ttR 2 c. Using R ~ 10 14 7?i4 cm, 3 times the Schwarzschild radius 
of a black hole of mass Mbh = 10 8 Mq, the magnetic held is B = 10 ^ (rjxL^i) 1 / 2 Rfl G. As¬ 
suming protons are accelerated in the same region, their maximum energy can be E max = 
eBR = 2.4 x 10 19 (r7x^44) 1 ^ 2 eV. This is problematic for the X-ray luminosities of Seyfert 
galaxies in Table 5, which are correlated with > 100 EeV cosmic rays, and the magnetic 
energy density must exceed the nonthermal X-ray energy density by a factor > 10 for pro¬ 
ton acceleration to ~ 10 2 ° eV. This additional energy could be accommodated if a sizable 
fraction of the Eddington luminosity, L E dd = 1-3 x 10 16 (Mbh/10 8 MQ) erg s _1 , could be 
converted to magnetic energy. The required cosmic-ray luminosities (k = 2.1) in Table 5 for 
Seyfert galaxies are above the Eddington luminosity for 2MASX J16311554+2352577, 2MASX 
J19471938+4449425, 2MASX J23272195+1524375. In case of Cygnus A, M hh and U Edd are 
an order of magnitude larger. The k = 2.3 cosmic-ray luminosities are more problematic. The 
opacity for photomeson (py) interactions with ex = 1 keV X-ray photons and the subsequent 
> 15 TeV neutrino production opacity is r P7 ~ 1 L^Rfl{ex/^ keV) -1 . 
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A hint of correlation that we found between the IceCube cosmic neutrino events and 
UHECRs with energy > 100 EeV should be investigated further by the experimental collabo¬ 
rations. Establishing a concrete correlation will be a ground-breaking discovery. Also a future 
extension of IceCube to increase its sensitivity in the > 1 PeV range will be very useful to 
probe the cosmic neutrino spectrum and if there is a cutoff in the spectrum. A cutoff in the 
spectrum is not natural at the PeV scale if the same sources produce > 100 EeV cosmic rays 
and neutrinos. Whether the weak AGNs, which are plentiful in the nearby universe, are the 
sources of UHECRs and neutrinos or not is a question that will continued to be debated and 
investigated in the years to come. 
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